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We study the pseudo-scalar, vector and axial current correlation functions in SU(2)-NJL model 
with scalar and vector couplings. The correlation functions are evaluated in leading order in number 
of colors N c . As it is expected in the pseudoscalar channel pions appear as Goldstone bosons, and 
after fixing the cutoff to reproduce the physical pion decay constant, we obtain well-known current- 
algebra results. For the vector and axial channels we use essentially that at spacelike momenta 
the correlation functions can be related to the experimentally known spectral density via dispersion 
, relations. We show that the latter imposes strong bounds on the strength of the vector coupling 

■ in the model. We find that the commonly used on-shell treatment of the vector and axial mesons 
ON \ (identified as poles at large timelike momenta) fails to reproduce the behavior of the corresponding 

correlation functions at small spacelike momenta extracted from the physical spectral density. The 
parameters of the NJL model fixed by the correlation functions at small spacelike momenta differ 
noticeably from those of the on-shell treatment. 
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The correlation functions of different colorless hadronic currents built of quarks and gluons offer a general tool 
. to study the structure of the QCD vacuum and the hadron spectrum Q. One can obtain information about the 
y—i 1 physical intermediate states from the pole structure of the correlation functions in the timelike region as well as from 
its behavior at spacelike momenta. For an exact theory both ways should lead to identical results. This, however, is 
not necessarily true for the effective models, which in the present lack of non-perturbative solution to QCD attract an 
increasing interest. Among them, the Nambu - Jona-Lasinio model, proposed (2| in early sixties in analogy with the 
r"| ■ BCS theory of superconductivity, plays an essential role. This model incorporates the basic symmetries of QCD and 
O i in particular, the chiral symmetry. It is able to manifest the spontaneous chiral symmetry breaking and, as a natural 
JL , consequence, the appearance of the pseudoscalar mesons as Goldstone particles. As such it offers a simple but quite 
successful working scheme to study the role of the dynamical chiral symmetry breaking mechanism in hadron physics 
(for review see refs. ^-Ql). However, the model suffers from two important drawbacks. First, it is not renormalizable 
L* ' and to make the theory finite one needs a finite cut off fixed in a consisitent way (usually to reproduce the pion 
decay constant) and the limit of large N c (number of colors). The typical value of the cut off is about 1 GeV. The 
second problem is the lack of confinement. Obviously, these drawbacks make the so far used on-shell treatment of the 
(axial) vector mesons questionable [§-[l0 12 In the present work we will show that the proper way to overcome 



this problem is to consider the corresponding current correlation function at small spacelike momenta. In this case 
we do not face the problem of the strong coupling to the qq continuum and we apply the model at small momentum 
region where it is designed for. To this end we make essentially use of the fact that the correlator at small spacelike 
momenta can be related via dispersion relations to well-established rich phenomenology without detailed knowledge 
of the dynamics at large timelike momenta. We will show that such a scheme allows to describe the vector modes in 
the meson sector of the NJL model in a consistent way. 

The simplest SU(2) Nambu- Jona-Lasinio (NJL) model is defined by lagrangean 

Cnjl = * (ip- m ) 9 + ^ [(**) 2 + (^r 75 *) 2 ] 

- ^ [{^l^ f + (<Ff 7M75 vI/) 2 ] , (1) 
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which contains chirally invariant local four-fermion interactions with scalar G s and vector G v couplings^. The current 
quark mass mo stands for both up and down quarks (assumed to be degenerated in mass). 

Applying the well-known bosonization procedure |14[ ] we arrive at the model generating functional expressed as a 
functional integral over auxiliary boson fields S, P, V^, A^ : 

f - - JVoTrlogf D[S,P,V U ,X U ] )+I m [S,P,V u ,A u ] , 

Z NJL = j VSVPVV^VA^e S V M V ' ? (2) 

where the quarks are already integrated out. To simplify our notations we adopt for P, and that P = P ■ f 
etc. The euclidean Dirac operator D is given by 

D = -id^ +S + iP>y 5 + + ,4 m 757m (3) 
and the trace over colors is explicitly done. The symmetry breaking term is included in the boson part: 

Im = I d M 2^i ( * 2 + p2) + *k ( ^ 2 + ^ - ■ (4) 

In the large N c limit the integral in eq.(||) is given by its saddle point value. For the translational invariant case 
(vacuum) we have a non-trivial stationary meson configuration given by 

S = M, P = 0, V„=0, 1^=0, (5) 

which means that the chiral symmetry is spontaneously broken. The quarks acquire a constituent mass M which is 
related to the scalar coupling constant G s via the stationary condition (saddle point): 

M = 2G S MN C Tr(-^gj) + m = m - G s < > . (6) 

Here < tyty > is the chiral quark condensate. As it was mentioned the model is not renormalizable and in order 
to make the fcrmion determinant finite it must be regularized. Thus, the model has four parameters, namely two 
coupling constants G s and G v , the current mass mo and the cut off A. Usually they are fixed reproducing the physical 
pion decay constant and the physical masses of pion and rho meson. It leaves one parameter free which is commonly 
chosen to be the constituent mass M because with eq.(||) one can eliminate G s in favor of M. In principle, one can 
use the empirical value for the quark condensate < > to constrain M. However, being quadratically divergent 
< ^*i> > is very sensitive to the details of the regularization scheme. 

In this paper we suggest an alternative way of fixing the model parameters G s , G v , m and A. To this end we 
study the correlation functions of the pseudoscalar, vector and axial-vector currents: 

P a {x) = *(a0i75y*, (7) 
V«{x) = *{x) 1 ^, (8) 

AJ(s)=*(s)7^y IS y*, (9) 

in the NJL model. The correlation functions arc defined by the vacuum-to-vacuum matrix elements of the currents 
and can be expressed in terms of invariant functions as follows: 

U P (Q 2 )= [ d i xe tqx {0\T{P a (x)P a (0)}\0) , (10) 



1 This lagrangean with vector couplings was considered first by Kikkawa 



2 



- q 2 9^v(Q 2 ) = J d 4 xe^(0 \T{V^(x)V;(0)}\ 0) , (11) 

(g„g* - ?V)nl(g 2 ) + g^X(Q 2 ) = | rW«*(o |T{^(x)^(o)}| o) , (12) 

where as usual Q is defined as Q 2 = — q 2 . In contrast to the vector current in the real world the axial- vector current 
is not conserved and the correlation function in the axial channel has an additional longitudin al p art 11^ (Q 2 ) which 
vanishes in the chiral limit too — > 0. In the NJL model[] the matrix element of rhs of eqs.f liofjl^ ) can be expressed 
as a path integral using the model generating functional (^). In the leading order in N c only two diagrams (shown 
in fig. |l|) contribute to the matrix elements. The first one is a simple one-quark loop whereas the second includes a 
boson line. For the pseudoscalar function we get the following result: 



%(Q 2 )/ a (Q 2 



<jng> 



n a (Q 2 ) = - - M - , (13) 

Gsfa(Q 2 )Z p (Q 2 ) Q 2 + JTgI Bp (Q*)f a {Q'>) 



where 



/a(Q2) 1 + 4M 2 G V Z P (Q 2 )- (W) 

Function Z p (Q 2 ) corresponds to the quark- loop diagram with two pseudoscalar-isovector vertices (»75T a ). As already 
mentioned, the model is not renormalizable and in order to make the quark loop finite we need a regularization. In 
our calculations we use the proper-time as well as the Pauli-Villars regularization scheme. Both schemes preserve the 
symmetries of the model. Here we present only the proper-time regularized expression for Z p (Q 2 ) Jl2[: 



+ 1 oo 
-1 A- 2 



M 2 + ^(l-« 2 ) 



(15) 



and for the condensate: 



N C , , f ds _ sM . 



<**>=-^Af / — e~ sM . (16) 

Details as well as regularized expressions for the quark condensate and the function Z p (Q 2 ) in the case of Pauli-Villars 
regularization can be found in ref . []l5| . 

In the chiral limit the pseudoscalar correlation function (O) develops the expected Goldstone pole at Q 2 — 0: 



n a (Q 2 ^o)^ 02 , 2 (17) 



< \E"I> > 5 

where the pion decay constant / 2 is given by the residue of the pole: 

f 2 = M 2 f a (0)Z p (0). (18) 

The non-zero current mass mo shifts the pion pole in eq.(|l7|) to Q 2 — —m 2 and one recovers the Gell-Mann - Oakes 
- Renner relation: 

/ 2 m 2 = -to < > +0(m 2 ) . (19) 
In the vector channel the correlation function H V (Q 2 ) is given by 

^■l iA 1 (20) 



2 Using different techniques the invariant functions are evaluated also in ref. JTJ| 
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where, similarly to Z p (Q 2 ), the function Z V (Q 2 ) corresponds to the quark-loop diagram with two vector vertices 
(7 M T a ) (lhs of fig. |l|). In particular, in the proper-time scheme it has the following form: 



M 2 + %(l-u 2 ) 



+ 1 oo 
-1 A- 2 

For the transverse invariant function of the axial current correlator in the NJL model we have: 



(21) 



2 TT' (G 2 \ = 1 Z v {Q 2 )Q 2 +±M 2 Z p (Q 2 ) 

W ,W ) 4 l + G^(Q 2 )Q 2 +4G„Af 2 Z p (Q 2 ) ' K ' 



whereas the longitudinal one includes only a pion contribution: 



q 2 K(q 2 ) = - - f ^: lf; — ■ (23) 



The correlation functions (p0|) ,(22) possess (at some model parameter values) poles in the timelike region Q 2 < 
which is in fact the pole of the propagator of the auxiliary (axial) vector field. This is used in the so-called on-shell 
treatment §-|l| where one fixes the vector coupling constant G v adjusting the position of the pole to the physical 
rho mass. The other model parameters, namely the cutoff A and the current mass mo are fixed reproducing the 
physical pion mass and the physical pion decay constant from eqs. (|l8|) , (|l9|) , whereas the constituent quark mass is 
treated as a free parameter as usual. To illustrate the on-shell scheme we present in fig.^Jthe corresponding parameter 
values (index "pole" ) as a function of M for the case of the the Pauli-Villars regularization. Since the results with the 
proper-time regularization are almost identical with those of the proper-time ones we do not discuss them separately. 
As can be seen the behavior of the coupling constants and the cutoff as function of the constituent mass is not smooth 
- the value M = m^ cp /2 is a singular point where all curves show a kink. At M > m exp /2 the vector current 
correlation function has a real pole at Q 2 — —m 2 which is an indication for a bound state in the spectrum. However, 
because of lack of confinement a second peak in the qq continuum appears, which is not taken into account in the 
definition of the physical p field. Obviously, it is hard to relate the second peak to the excited p(1450)-resonance. 
At smaller values of M < m exp /2 in the vector channel a broad resonance appears in the qq continuum because the 
corresponding pole moves to the complex plane of Q 2 . The position of the peak depends on the vector coupling 
constant G v . Using thisQ we determine G v fixing the position of the peak in the vector spectral function ImTL at 
the physical p mass. However, as can be seen from fig.|| it is not always possible to find a solution for G v : for M 
between 230 and 260 MeV in fact there is no such a solution. It should be also noted that at M > m e p xp /2 both G v 
and G s , as well as the cutoff A show a stronger dependence on M as G v is dominant (G v > G s ). This changes for 
smaller values of M where G s and A stay almost constant whereas G v decreases. In the axial channel for the values 
of the constituent mass M considered (200 - 450 MeV) the A\ meson appears as a very broad peak being centered 
around 1 GeV in the qq continuum. In Fig.^j we also show the parameter values (index "grad" ) obtained by means of 
a gradient or heat-kernel expansion of the effective action | p^p^7| . This procedure is frequently used to fix the vector 
coupling constant G v in NJL model and in fact it is an approximation to the on-shell description of vector modes. 
As can be seen from Fig.^ this procedure provides a quite crude approximation to the on-shell values because of the 
large masses of the vector mesons. From the above discussion one concludes that the on-shell treatment is not able 
to provide a consistent description of the vector meson modes in NJL model. 

As a next step we suggest to make use of the fact that at spacelike momenta Q 2 > the vector and axial-vector 
current correlation functions can be related to the experimentally known spectral density via dispersion relations. In 
the vector channel the invariant function H V (Q 2 ) satisfies a dispersion relation with one subtraction: 

n.^-n.m-f/d.^^, (24) 

where the physical spectral density Imll„(s) is experimentally accessible: 

ImlLj(s) = — . (25) 



3 Such a procedure was used in refs. 
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For the rhs we use a simple (p + continuum) parameterization 1T3] of the experimental data 

Km 2 n 1 / a = \ 

Imn„( S ) = —£6(8 ~m 2 p ) + —(l + ^ 8 s - s , (26) 



with the parameter values: 



2 

= 0.77 GeV, (^=2.36, s = 1.5GeV 2 . (27) 



Eqs.(p4[) and (^6|) provide an independent way to fix the vector coupling constant G v confronting the model with the 
experiment by equating the phenomenological dispersion integral and the model correlation function at low positive 
Q 2 : 



n -(°)- v / ds 



lmn v (s) _ 1 Z V (Q 2 ) 



5 (s + Q 2 ) 41 + G v Z v {Q 2 )Qi 



The resulting G v as well as A and G s as a function of M, obtained by adjusting (x 2 -fit) the lhs of eq. ( p8[ ) to the rhs 
with the physical spectral density ( p6| ) at low positive momenta (we consider < Q 2 < 0.5 GeV 2 ), are shown (index 
"corr") in fig.^j a) and b). Apparently the correlation function at spacelike momenta and the on-shell treatment yield 
rather different results. In contrast to the curves resulting from the on-shell scheme those of the correlation function 
have the advantage to behave smoothly as a function of M. This, and the very fact that the NJL model is designed 
for low-energy structures involving small momenta, makes the above correlation method preferable. With increasing 
M larger values for G v are needed to fit the phenomenological side whereas the G s is slightly decreasing. 

In the axial channel, since the longitudinal part of the axial correlation function is almost independent on the G v 
and M, we will concentrate our further discussion on its transverse part 11^. It obeys a dispersion relation with two 
subtractions: 

Similar to the vector channel for Imll^Q 2 ) we take a simplified (A\ + continuum) parameterization 

1 / CK 



7T 



9 P 1 8ir\ - ' 



m\ — 1.26 GeV, ( -^-) = 0.15 - 0.18, s = 1.7GeV 2 . (31) 



with parameters llf 



Confronting the model Ai-correlation function fl22|) and the phenomenological expressions (|29|), (|30| ) one has another 
way to fix the vector coupling constant G v . The obtained vector coupling constant (G^F^ ) as a function of M is shown 
in figj^ a) in comparison with the results from fitting the p phenomenology as well as from the on-shell p treatment and 
the derivative expansion. In contrast to the vector channel the model axial function 11^ (Q 2 ) depends much stronger 
on the constituent mass and because of that the model is able to fit simultaneously both p- and ^i-channels only 
within a narrow window of values for the constituent quark mass M around 240 MeV. It is interesting to note that 
in this window the on-shell values are also very close to those obtained from the correlation functions. The particular 
value of M depends on the parameter values (27),(|3l"|) used for the (resonance + continuu m) par ameterizations. In 



fact the estimate of the constituent mass M « 300 MeV in the random instanton model |21jJ22| of QCD vacuum, 
applied successfully also to the meson correlation functions |23] ], is not far from our numbers. The resulting vector 
coupling constant G v is of the same order as the scalar one G s . In this point we disagree with ref. |20[ in which 
constraining the NJL model via QCD sum rules it is concluded that G v should be almost zero or at least G v <G S . 
It is due to the fact that in contrast to us in ref. [^0| the NJL model is treated in a quite simplified approximate way 
ignoring the full Q 2 -dependence of the correlation functions as well as the relations between the model parameters 
fixed to reproduce the physical pion properties. 

In our calculations so far we used simple parameterizations (p6|) , (|30|) for the experimental spectral densities in 
which the resonance widths are neglected. In order to check this approximation we repeated the calculations using 
a much more elaborated (finite-width resonance+continuum) parameterization plf of the experimental data. The 
obtained results follow qualitatively the picture of fig. 0. The only difference is that the narrow window, where in 
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both vector and axial channels the experimental low-energy behavior of the correlation functions can be reproduced 
simultaneously, is shifted to higher values of the constituent mass M w 260 MeV. 

We also related the model with G v , fixed in the present scheme, to the chiral effective lagrangean of Gasser and 
Leutwyler ^6|. Following ref. [^5| we calculated the corresponding low-energy coefficients l\ — 1$ in the notation of 
ref. |p6| and low-energy pionic characteristics. As can be seen from Table | our theoretical predictions are in good 
agreement with the empirical values. 

To summarize, in this letter we study the vector meson modes in the NJL model in terms of current correlation 
functions. We find that the on-shell treatment of the vector meson fails to reproduce the low-energy behavior of 
the correlation functions fixed by the experimentally known spectral density via dispersion relations. We offer a 
different scheme to fix the NJL-model parameters based on the behavior of the correlation functions at small spacelike 
momenta. Treating the vector coupling constant as a free parameter the model is able to reproduce the phenomenology 
in both the axial- and vector channels for the constituent mass only in a narrow window around 240 MeV. The vector 
coupling constant is by no means zero and is of the same order as the scalar coupling constant. In principle, these 
two schemes, the on-shell treatment (large timelike momenta) and the one based on the behavior of the correlation 
function at small spacelike momenta, should provide identical results within an "exact" theory, which however is not 
the case for an effective model like NJL. In the latter case the correlation function method appears to be preferable. 
One also should keep in mind that the present considerations are done in leading order in N c . The inclusion of 1/N C 
quantum boson (loop) corrections could in principle change the present results. 
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FIG. 1. Diagrams contributing to the correlation functions in leading N c order. 

FIG. 2. The parameters of the model as function of the constituent mass M: a) vector coupling constant G v fixed in 
the on-shell treatment G% ole , from derivative expansion G 9 v rad , and from the correlation functions in the case of p- and 
ai-phenomenology fits; b) the same as a) but for the scalar coupling constant G s and the cut off A. 



7 



TABLE I. Low-energy SU(2) coefficients and low-energy pion properties for constituent mass M = 235 MeV. 



SU(2) coefficients 


model 


ph.en.orn. 


TV properties 


model 


exp. 


A [MeV] 


1400 




ail 


0.16 


26 + 05 


— <r rrn ^ 1 / 3 [MpVl 
\ yy ^ [ivic v j 


OIU 


ZiXJO ZL H:\J 


°o 


n 1Q 


n oc _i_ o no 


mo [MeV] 


3.5 


7 ± 3 


„2 
da 


-0.045 


-0.028 ± 0.012 


h 


-4.17 


-4.56 ±1.5 


hi 


-0.084 


-0.082 ± 0.008 




3.02 


2.3 ±0.7 


a\ 


0.035 


0.038 ± 0.002 


h 


2.45 


-1.0 ±3 


a% x 10 4 


9.8 


17.0 ±3.0 


U 


2.44 


1.1 ±0.4 


a\ x 10 4 


-1.1 


1.0 ±3.0 


h 


9.1 


10.4 ± 1.3 


< r 2 >J [fm 2 ] 


0.43 


0.52 ±0.15 


h 


12.7 


13.1 ± 1.3 


< r 2 >Z [fm 2 ] 


0.37 


0.44 ± 0.03 
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